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ABSTRACT 


It  is  noted  that  certain  common  linear  wave  operators  have  the 
property  that  linear  variation  of  the  initial  data  gives  rise  to  one- 
dimensional  evolution  in  a  plane  defined  by  time  and  some  direction 
in  space.  The  analysis  is  given  for  operators  arising  in  acoustics, 
electromagnetics,  elastodynamics,  and  an  abstract  system. 
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1  Introduction 


In  this  paper  we  point  out  an  interesting  fact  concerning  three  common  instances  of  linear  wave 
propagation  (acoustic,  electromagnetic,  and  elastic).  In  each  case,  we  consider  the  initial-value 
problem  in  some  region  of  space  sufficiently  small  that  the  gradients  in  the  data  can  be  taken 
as  constant.  Then  the  evolution  of  that  data  in  time  corresponds  precisely  to  that  produced  by 
one-dimensional  waves  travelling  in  some  distinguished  direction.  In  the  three  instances,  this 
direction  is  defined  *  by  the  pressure  gradient,  by  (V  x  5)  x  (V  x  D)  (not  the  Poynting  vector), 
and  by  any  vector  perpendicular  to  V  x  u.  We  also  consider  an  abstract  system,  to  which  all 
2x2  systems  in  a  certain  class  are  equivalent. 

2  Statemert  of  Result 

Consider  a  set  of  linear  partial  differential  equations  in  n  unknowns  and  d  space  dimensions, 

i=d 

5<q+5Z>f.5*riq  =  0.  (1) 

t=i 

A  simple  wave  solution  of  (1)  is  a  solution  of  the  special  form 

q(£,t)  =  Tif{C  x-  At) 

where  ^  is  a  unit  vector  in  the  direction  of  wave  propagation,  A  is  the  speed  of  propagation,  and 
/  is  an  arbitrary  scalar  function.  Such  solutions  exist  if  and  only  if  tt  is  a  right  eigenvector  of 

i^d 

Ae  =  '£eiAi. 

«=i 

Linearly  varying  initial  data  for  (1)  can  be  considered  as  a  vector  space  V  in  R"*^,  conveniently 
displayed  as  d  column  vectors 

D  —  [dxiq>  dxjq*  ‘  *  ■] 

Data  giving  rise  to  a  simple  wave  is  of  the  special  form 


Di  = 

Denote  by  St  the  subspace  of  V  containing  such  data,  with  I  fixed  and  A  ^  0,  i.e.,  the  set  of 
data  giving  rise  to  one-dimensional  wave  motion  in  the  direction  t. 

The  evolution  operator  A  that  produces  dtq  from  an  element  of  P  is  an  n  x  nd  matrix  whose 
nuUspace  A/’(A)  comprises  the  set  of  linearly  varying  steady  solutions  to  (1).  We  claim  that  for 
some  interesting  evolution  operators,  such  as  those  mentioned  in  the  Introduction,  an  arbitrary 
element  of  T>  can  be  represented  as  a  sum  of  elements  lying  in  Af  and  elements  lying  in  some  St, 

D  =  J2PjNj  +  J2<^kSt,k  (2) 

i  * 

‘The  notation  used  in  this  paper  is  that  bold  type  represents  a  vector  such  as  a  set  of  unknowns;  vectors  with 
geometric  meaning  in  have  arrows  surmounting  regular  type 


Note  that  we  do  not  claim  that  Af  +  St  =  T>,  but  that  (when  the  trick  works)  the  direction  / 
can  be  chosen  so  that  Af  +  St  includes  any  given  element  of  V.  To  evaluate  the  coefficients  in 
(2)  we  operate  on  both  sides  with  >1,  giving 

=  52  ^k^krt,k  (3) 

k 

which  is  merely  a  decomposition  of  the  time  derivative  onto  the  eigenvectors  for  wave  motion 
along  1.  This  differs  from  the  one-dimensional  decomposition  (Riemann  problem)  in  that  the 
direction  of  t  is  unknown  and  enters  into  the  equations  nonlinearly.  In  the  following  sections  we 
see  how  this  works  out  for  various  special  cases. 


3  Examples 

3.1  Examples  with  Unique  Solutions 

3.1.1  Acoustics 

For  the  equations  governing  acoustic  waves  in  a  uniform  medium  we  have 

0  It  t2  ^3 

0  0  0 
"  4  0  0  0 

4  0  0  0 

where  the  unknowns  are  the  pressure  and  velocities,  and  units  have  been  chosen  to 

make  the  sound  speed  unity.  The  eigenvalues  are  Ai  ,2,3,4  =  -1,0,0, 1  and  the  non-stationary 
eigenvectors  are 

ri  = 

The  expansion  (3)  is 


which  can  only  be  solved  if 
when  the  wavestrengths  are 


t  \\  vt  =  -gradp. 


1  -  1- 
at  = -{t  ■  vt  -  Pt),  Q4  = '  vt  +  pt). 


(5) 


3.1.2  Electromagnetics 

Maxwells  equations  in  free  space,  for  the  magnetic  field  B,  and  the  electric  field  D,  in  units  for 
which  the  speed  of  light  is  unity,  give  rise  to 
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The  eigenvalues  Ai, 2, 3, 4, 5,6  are  -1, -1,0, 0,1,1  and  the  non-stationary  eigenvectors  are 

s  t  t 

ri  =  -  ,  r2  =  ^  ,  rs  =  «  ,  re  = 

t  s  —s 

where  s,  t  are  any  pair  of  unit  vectors  such  that  t,  s,  t  form  a  right-handed  orthogonal  system. 
The  expansion  (3)  is 

=  ai  -  +02  -  +0(5  -  +0(6  r  •  (') 

X  S  S  i 

implying 

Bt  =  (ae  -  Q(i)s  +  (05  -  02)^ 

Dt  =  -(ae  +  0(1  )t  —  (0:5  +  02)5 

For  a  solution  to  exist,  both  Bt  and  Dt  must  lie  in  the  plane  spanned  by  s  and  t,  which  is  the 
plane  normal  to  i.  Then 

£*  II  Bt  y.  Dt  =  curlB  x  curlD  (8) 

and  we  easily  find 

0(1  =  -BfS-Dft, 

02  =  -Bft-  DfS, 

0(5  =  Bf  t-  Dt  •  s,  (9) 

oo  =  Bf  s  —  Dt  - 1. 

3.2  Examples  with  Nonunique  Solutions 
3.2.1  ‘Cauchy-Riemann’  Equations 

Consider  the  system  obtained  by  adding  time-dependent  terms  to  the  Cauchy-Riemann  equa¬ 
tions; 

«t  +  «3r  +  Wy  =  0 

Vt-Vx  +  Uy  =  0  (10) 

for  which  the  matrix  At  is 

(11) 

£2  — £1 

Gilquin,  Laurens  and  Rosier  [1]  have  shown  that  any  strictly  hyperbolic  2x2  system  can  be 
reduced  to  this  case  by  a  transformation  of  variables,  and  Noelle  [2]  finds  that  any  linear  initial 
data  is  the  sum  of  three  simple  waves. 

The  eigenvalues  of  j4/  are  Ai  =  -1,A2  =  1,  with  eigenvectors 


Bt 

Dt 
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These  can  be  normalised  and  rewritten  in  terms  of  the  wave  direction  B  as 


ri  = 

and  the  expansion  (3)  becomes 


—  sin  0/2 

cos  0/2 

cos  0/2 

»  *'3  = 

sin  0/2 

«< 

vt 


=  oi 


sin  6/2 
—  cos  6/2 


+  Otj 


cos  6/2 
sin  0/2 


(12) 


This  can  be  solved  for  any  choice  of  6:- 

oi  =  ut  sin  6/2  —  Vt  cos  0/2,  aj  =  «<  cos  0/2  +  sin  0/2.  (13) 

3.2.2  Elastodynamics 

The  equations  governing  elastic  waves  in  a  uniform  isotropic  medium  are  of  the  form  (1)  with 
q  =  (n,p,w)^  where  v  is  the  velocity,  p  is  the  trace  of  the  strain  tensor,  and  w  forms  the 
antisymmetric  (rotational)  part  of  the  strain  tensor.  The  strains  are  related  to  the  displacements 
u  by 

p  =  div  «, 


and 


u  =  curl  tt. 


The  matrix  Ai  can  be  written  as 
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(14) 


Here  units  have  been  chosen  so  that  transverse  waves  (S-waves)  travel  with  unit  speed,  and 
longitudinal  waves  (P- waves)  with  speed  \/b,  where  the  parameter  b  is 

I.  ,  2M 

and  L,  M  are  the  Lame  coefRcients. 

The  matrix  (14)  has  eigenvalues  Ai, 2, 3, 4, 5, 6, 7  =  — >/&,-!  -  1,0,  l,l,\/5,  with  non-stationary 
eigenvectors 


Vb( 

S 

ri  = 

-1 

0 

,  ^2  = 

0 

-t 

,  tz  = 

t 

3 

rs  = 

0 

,  re  = 

0 

,  rr  = 

3 

t 

t 

0 

—s 

y/hi 

1 


(15) 

(16) 
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Thus  the  expansion  (3)  for  this  case  is 


Vt 

-U 

Pt 

= 

Vb 

+  0(2 

Qt 

0 

-t 

t 

+  0(3 

0 

+  as 

0 

s 

5 

S 

u 

+  0(6 

0 

+  07 

yfh 

t 

0 

(17) 


From  this  we  observe  that  Qt  lies  in  the  plane  spanned  by  s,  t,  so  that  I  must  be  chosen  normal 
to  wt  =  curltT.  For  any  such  choice  of 


a?  +  0(1 

= 

PjlVb, 

07  -  ai 

= 

{•  vt/b, 

0(6  +  0(2 

0(6  —  0(2 

= 

S-Vt, 

0(5  +  O3 

= 

£-wi. 

Os  -  0(3 

= 

t  •  Vj. 

(18) 


3.3  An  Example  with  no  Solution 

The  following  is  based  on  a  suggestion  of  Jonathan  Goodman  (Courant  Institute,  NYU).  Con¬ 
sider  a  system  comprising  two  copies  of  the  acoustic  equations,  with  unknowns  (wi,pi,U2,p2). 
The  resulting  one-dimensional  wave  problem  has  a  solution  only  for  data  with  ||  V2,f 


4  Applications 

These  results  have  potential  application  to  the  numerical  analysis  of  wave  motion,  where  the 
solution  is  often  assumed  to  be  linear  within  elements,  and  where  the  update  strategy  may 
depend  on  whether  waves  are  entering  or  leaving  some  region.  For  a  critique  of  methods  where 
wave  directions  are  determined  relative  to  ceU  boundaries,  see  [3].  Dissatisfaction  with  such  an 
approach  leads  to  an  interest  in  representing  multidimensional  waves  in  a  way  that  is  coordinate- 
free,  but  simple  enough  for  computation.  The  lack  of  a  general  pattern  may  not  be  a  drawback; 
exploiting  behaviour  that  is  very  specific  to  a  particular  equation  set  may  be  the  most  effective 
approach. 

Specific  applications  will  turn  on  details  of  numerical  technique  as  well  as  the  analytical 
decompositions;  for  that  reason  they  are  not  explored  here,  but  reference  may  be  made  to 
computations  of  the  compresible  Euler  equations  by  splitting  finite-element  gradients  into  steady 
and  unsteady  components  [4],  and  to  a  boundary  condition  for  Maxwell’s  equations  that  uses 
the  above  formulae  to  identify  the  direction  in  which  waves  exit  the  boundary  [5]. 

In  a  non-numerical  context,  there  is  precedent  for  this  analysis  in  the  work  of  Frohn  [6]  who 
studied  steady  three-dimensional  supersonic  flow  by  the  method  of  strained  coordinates,  where 
the  objective  is  to  find  a  coordinate  perturbation  such  that  the  singularities  of  linear  theory 
become  realistic  shocks  and  rarefactions.  She  selected  a  plane  containing  the  streamline  and  the 
pressure  gradient  as  the  one  in  which  a  two-dimensional  straining  should  be  applied 


5  Conclusions 

Given  a  set  of  hyperbolic  partial  differential  equations  in  D  space  dimensions  (D  >  1)  and  time, 
together  with  linearly  varying  initial  data,  one  may  seek  directions  in  which  the  evolution  of  the 
problem  takes  place  one-dimensionally.  We  have  given  examples  in  which  the  search  succeeds 
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(and  a  rather  artificial  counterexample).  Since  a  practical  computational  algorithm  can  only  be 
based  on  a  finite  number  of  wave  directions,  this  offers  a  possible  basis  for  the  creation  of  such 
algorithms. 
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